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Abstract. The processes contributing to the evolution of 
an initially weak magnetic field in a differentially rotat- 
ing star arc reviewed. These include rotational smoothing 
(akin to convective expulsion) and a list of about 5 insta- 
bilities, among them magnetorotational instability, buoy- 
ancy instability, and pinch-type instabilities. The impor- 
tant effects of thermal and magnetic diffusion on these 
instabilities are analyzed in some detail. The first instabil- 
ity to set in is a pinch-type instability. It becomes impor- 
tant in modifying the field configuration before magnetic 
buoyancy-driven instabilities set in. The evolution of an 
initially strong field remains a more open question, includ- 
ing the old problem whether dynamically stable magnetic 
equilibria exist in stars. 
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1. Introduction 

A number of different processes cause stars to rotate differ- 
entially. Stars with convective envelopes experience brak- 
ing of their rotation by a magnetic stellar wind (Schatz- 
man 1962). This torque acts on the convective envelope, 
and slows it down relative to the radiative core. Secondly, 
in an evolving star the core contracts and spins up, while 
the envelope expands and slows down. If core and enve- 
lope were to conserve their angular momentum separately, 
the core of a giant would end up rotating some 10 5 times 
faster than its envelope. Finally, the Eddington-Sweet cir- 
culation in a rotating star causes differential rotation by 
transport of angular momentum. In a steady state, such 
that dissipation and driving of the differential rotation by 
this circulation balance each other, the rotation decreases 
outward by some 30% between the core and the surface 
(Zahn 1992). 

Opposing this differential rotation are friction pro- 
cesses due to hydrodynamical instabilities or magnetohy- 
drodynamical processes. The strength of such friction de- 
termines how much rotation can be left in the end prod- 



ucts of stellar evolution, white dwarfs and neutron stars. 
If they are effective enough to maintain an approximately 
uniform rotation in giants, for example, the observed rota- 
tion of white dwarfs and pulsars is not a leftover of the ro- 
tation of their progenitors (Spruit & Phinney 1998, Spruit 
1998). Explosion mechanisms of supernovae that rely on 
rotation (Bisnovatyi-Kogan 1970, LeBlanc & Wilson 1970, 
Meier et al. 1976, Rampp et al. 1998), and explanations 
of the morphology of planetary nebulae and objects like 
r\ Carinae in terms of stellar rotation (Heger & Langer 
1998) also depend crucially on the degree of differential 
rotation that can persist in a stellar interior. The rele- 
vant hydrodynamical processes have been classified (Zahn 
1974, 1983). Their net effect on angular momentum trans- 
port is somewhat uncertain and in stellar evolution calcu- 
lations is usually parametrized with arbitrary adjustable 
coefficients. 

An important test case for such parametrizations is 
the internal rotation of the Sun, which has been measured 
with remarkable precision though helioseismology (Cor- 
bard et al. 1997, Schou et al. 1998). The most striking 
result of these measurements is the near-uniform rotation 
of the radiative core of the Sun. Over a depth range of only 
O.O5i?0, the differential rotation in the convective envelope 
(30% faster at the equator than at the poles) changes to a 
state of uniform rotation in the core. Though the measure- 
ments are less accurate in the inner core (t/Rq < 0.3) and 
near the rotation axis, no deviation from uniform rotation 
has yet been detected in the core. 

It has been known for a long time that such a low 
degree of differential rotation is incompatible with the 
currently known hydrodynamic transport mechanisms of 
angular momentum. Circulation is ineffective for a slow 
rotator like the Sun. The known hydrodynamic instabil- 
ities are also unlikely to contribute, since the conditions 
for their occurrence are not satisfied in a core rotating 
as slowly and uniformly as the present Sun (Spruit et al. 
1983). 

It is therefore natural to look for magnetohydrody- 
namic mechanisms of angular momentum transport in sta- 
bly stratified layers of stars. This may sound like a daunt- 
ing prospect, given the reputation of magnetic fields for 



2 



H.C. Spruit: Differential rotation and magnetic fields in stellar interiors 



complexity and lack of firm results. The situation in the 
context of the present problem is in fact quite promising. 
Over the past 4 decades, a substantial MHD literature 
bearing on the problem has developed. In this paper I first 
review the various contributing processes, conclude that 
pinch-type instability of an azimuthal field (Tayler 1973) 
is of particular relevance. Then I show that the effects of 
magnetic and thermal diffusion arc quite important, and 
how they can be quantified using Acheson's (1978) disper- 
sion relation. 

1.1. statement of the problem 

Suppose we start with a differentially rotating star with 
rotation rate f2(r) and field configuration Bo- We wish 
to know how the field and the rotation evolve in time, 
and to answer questions like: is the field able to make the 
star rotate uniformly? If so on what time scale? What 
kind of field configuration remains at large time? Though 
these questions can not be answered fully because of a few 
missing pieces of theory, it turns out that a fairly detailed 
picture can be drawn of the magnctohydrodynamics of 
an initially weak field in a differentially rotating star. In 
most of the analysis, I will assume that there is no external 
torque on the star, and that its internal structure is not 
evolving. These assumptions are made for convenience and 
clarity in delineating the processes involved; they can be 
easily relaxed. 

A number of different MHD processes are involved in 
the problem stated. I introduce them by starting with a 
number of unrealistic simplifying assumptions and then 
relaxing these. The evolution of the field and the rotation 
depends on the strength of the initial field Bo. In most 
of the following, I assume the initial field to be weak, in 
the sense that the Alfven travel time through the stars 
is much longer than the rotation period. In terms of the 
Alfven frequency ft a, 



fit 



D 



RiAnt)) 1 / 2 ' 



(1) 



where p is the mean density and B a mean field strength, 
the assumption is that fix <C fl. For the current Sun, this 
is satisfied if the field is small compared to a megaGauss, 
which in all likelihood is the case. For early type stars, with 
their larger rotation rates, the assumption is even more 
likely to be justified. It is less well justified, however, in 
the cores of giants, if these corotate with their envelopes. 

2. Winding-up of weak fields (no diffusion, v$ only) 

2.1. winding-up 

Start with ignoring magnetic diffusion, so that the induc- 
tion equation is 



and by assuming that the only motions to be considered 
are axisymmetric, purely azimuthal motions, 



^0 = rfl(r, 9). 



(3) 



Thus we are ignoring for the moment the motions in the 
(r, 9) plane due, for example, to instabilities. Also ignored 
are viscosity and thermal diffusion. Consider first the case 
of an axisymmetric magnetic field (aligned with the rota- 
tion axis). In the usual way, an axisymmetric field B can 
be written in terms of a stream function ip: 



B r 



-dgip, B e 



r sin6* 



d r ip. 



(4) 



This function is constant along field lines. Each value of tp 
labels a different magnetic surface (surface generated by 
rotating a field line around the axis). 

The induction equation then has the components 



d t B r — d t B s 



0. 



d t B d 



rsin6»B p • Vfi. 



(5) 



Under the assumptions made the poloidal field component 
Bp = (B r ,Be 7 0) does not change in time, and the equa- 
tion for B^ can be integrated in time: 



Ba, = Nn ■ B 



where 



N(r,9,t) =rsin0 J | Vfl\ 



(6) 



(7) 



is (modulo a factor 2tt), the number of 'differential turns', 
or the rotational displacement, and 



n = VO/|VO| 



(8) 



is a unit vector in the direction of the gradient of fl. The 
torque per unit surface area, 



t = r sin 9B r B ( p /An 



(9) 



increases linearly with the displacement N, hence the 
torque causes the rotation to execute a harmonic oscilla- 
tion. Under the assumptions made, the oscillation remains 
linear even though the displacement can be very large. 
This is because of the linear behavior of a pure Alfven 
wave at arbitrary amplitude. The oscillation is, in fact, an 
Alfven wave traveling along the magnetic surface gener- 
ated by the poloidal field lines with a given value of the 
stream function ip. Its period is the Alfven travel time 
/ ds/VAp along this field line, where Va p is the Alfven 
speed based on the poloidal field strength (Mestel 1953). 

Considering the evolution of the field on a given time 
scale to (for example the age of the star since its forma- 
tion), we can define a critical field strength such that the 
oscillation period t\ = R/va equals t: 



3iB = Vx(ux B), 



(2) B M = (4TT P ) 1/2 R/t , 



(10) 
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where p is a typical density in the star, and R its radius. 
For example, take for to « 10 17 s the age of the present 
Sun, and for p its mean density « 1. Then Bm ~ 2pG 
(Mestel 1953). That is, if the initial field strength is more 
than a /iGauss or so, the field eventually gets wound up 
to such an extent that Lorentz forces start affecting the 
internal rotation of the Sun. 

2.2. phase mixing 

If the initial field is large compared with Bm, not a very 
strong condition, the rotation on each magnetic surface 
oscillates around a mean value given by the total angu- 
lar momentum on this surface. Since Alfven waves do not 
couple across magnetic surfaces, neighboring surfaces os- 
cillate independently, have different periods, and the os- 
cillations on them increasingly get out of phase with each 
other. The distance between points with oscillation phase 
differing by n, say, decreases as 1/t. After a finite time, 
this distance becomes short enough that magnetic diffu- 
sion starts becoming important. Thus, it is necessary to 
include magnetic diffusion into the picture. 

3. Winding up and oscillation with magnetic dif- 
fusion (v<f, only) 

3.1. Rotational smoothing 

We now relax the asumption that the initial field is ax- 
isymmetric, and also include magnetic diffusion. The in- 
duction equation is then 



d t B = V x (u x B) - ??V 2 B, 



(11) 



where for simplicity I have assumed that the magnetic 
diffusivity is independent of position. This is a good ap- 
proximation, since the diffusion effects encountered here 
are effective only on small length scales. The poloidal field 
can be decomposed into two components 



Bp — B a + B n , 



(12) 



where B a is the azimuthal average of B, i.e. c^B a = 0, 
and B n the non-axisymmetric part, whose azimuthal av- 
erage (B n )0 vanishes. Let these components initially be of 
similar strength. 

If the initial field is weak, the field lines of the initial 
poloidal field B p0 are wound up tightly before the restor- 
ing Lorentz forces become effective. Let 



q = r\vn\/n 



(13) 



be a dimcnsionless measure of the local rate of differential 
rotation. If the initial poloidal field is smooth and vary- 
ing on a length scale of the order of the radius of the 
star, q ~ O(l). Since the induction equation is linear, the 
winding up process can be considered separately for the 
components B a and B n . Since the azimuthal average of 



B n vanishes, it has opposite polarities as a function of 
azimuth at any given (r, 9), and the azimuthal angle be- 
tween opposite polarities is tt or less. Points with field of 
opposite polarity at the same 9 but separated by a small 
distance I in the direction of the rotation gradient, will be 
brought to the same azimuthal angle after a time t given 
by: 



tClql rts nr. 



(14) 



The time scale on which these opposite polarities will can- 
cel each other by magnetic diffusion 77 is 



t d = l 2 / V . 



(15) 



This cancellation leads to a decay of the poloidal field, 
governed by the approximate equation 



<9 t ln_B n « l/t d = 7] 
which integrates to 

(t/tnf 



ftVlq 
\ nr 



B 



where the rotational smoothing time is 

1/3 



tn = 



3rV^ 
i]fl 2 q 2 



(16) 



(17) 



(18) 



Because of the steep dependence on t, the non- 
axisymmetric component of the poloidal field effectively 
disappears after only a few times in. 

This process has been studied extensively in the con- 
text of the kinematic evolution of magnetic fields in con- 
vection, and is called convective expulsion there (Zel- 
dovich 1956, Parker 1963, Weiss 1966). In the present 
context of differentially rotating stars, it was studied by 
Radler (1980, 1986). The remarkable effect of differen- 
tial rotation on a weak field is thus to expel the non- 
axisymmetric field components from the star on a finite 
time scale. For example, for an initial Sun rotating at 
fio = 3 1CP 5 (ten times the present-day rate), q s=s 1, with a 
diffusivity 77 w 10 3 , tn is only « lOOyr. If nothing else were 
to happen, the net result of differential rotation would be 
to make the field axisymmetric on a time scale which is 
very short compared to its life time. 

This, however, applies only if the initial field is suffi- 
ciently weak, so that the winding-up can be treated as a 
kinematic process. A field is weak in this context if mag- 
netic torques, which affect the differential rotation, do not 
become effective during the rotational smoothing process. 
This is the case if the initial Alfven travel time 



t AQ = r(4np)^ 2 /B 

satisfies 

*ao > tn, 



(19) 



(20) 
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which translates into a condition on the initial field 
strength: 



Bq < Bi, 
where 

Bi = r{4np) 



1/2 ( Q 



2„2\ 1/3 



3r 2 7r 2 



(21) 



(22) 



For the same initial Sun, B\ sa 30G. In terms of the cor- 
responding Alfven frequency: 



n 



3vr 2 



1/3 



V 
Qr 2 



1/3 



(23) 



If the poloidal field satisfies (pT|), the non-axisymmetric 
field component is smoothed out on the time scale (|l~8|), 
after which there remains an axisymmetric field (Radler, 
1986). For this to be the case, the initial Alfven travel 
time has to be long compared with the rotation period 
by a factor (fir 2 /^) 1 / 3 . If the initial field is larger than 
(p2|), rotational smoothing does not happen and instead 
the differential rotation is damped out by the process of 
phase mixing (all of this still under the assumption that 
i>0 is the only velocity component). 

3.2. Phase mixing 



Assume first that condition (21) is satisfied, so that ro- 
tational smoothing is effective. How does the field then 
evolve after being axisymmetrized? Since the smoothing 
process is linear, the axisymmetric and non-axisymmetric 
field components evolve independently. The symmetric 
component continues to be wrapped up after the non- 
axisymmetric component has been smoothed away. Even- 
tually, if time permits, the azimuthal field will become 
strong enough to oppose the wrapping. Phase- mixing then 
starts acting. Oscillations on neighboring magnetic sur- 
faces get out of phase, and the length scales in the Alfven 
wave field become increasingly smaller, until magnetic dif- 
fusion becomes effective and damps out the wave [Spruit 
1987, Roxburgh 1987 (private communication)]. If the 
Alfven travel time varies through the star by a factor 
qA = A£a/£a, the phase difference between magnetic sur- 
faces separated by a distance I is of the order flAtqAl/r. 
After a time t this phase difference reaches a value of it if 
the length scale I is 



I 



QA^At 



(24) 



The magnetic diffusion time on this length scale is = 
I 2 jr\. The wave amplitude A gets damped by magnetic 
diffusion on this time scale, 



d t lnA = l/t d = r\ 



V 7rr 



(25) 



This is the same equation as in the case of rotational 
smoothing, but with the wave frequency Qa replacing the 
rotation rate fl. The wave amplitude decays like 

A = A Q e- {t/t ^\ (26) 

where the phase mixing time scale t p is 

^■W • <27) 

When the oscillations have damped out, the rotation rate 
is constant on each magnetic surface: 



n = f(ip), (t>t p ). 



(28) 



The phase mixing process in Alfven waves also plays a role 
in other astrophysical situations, see Heyvaerts & Priest 
(1983), Sakurai & Granik (1984), Petkaki et al. (1998). It 
is analogous to the phase mixing process in plasma physics 
only in a rather broad sense. 

Expression ( p7| ) for the time scale is very similar to 
that for rotational smoothing (|l8|), but its consequences 
are rather different. Rotational smoothing removes the 
non-axisymmetric components of the magnetic field, phase 
mixing damps out the variations of Q along magnetic sur- 
faces. This has been studied with numerical examples by 
Mestel, Moss & Tayler (1988, 1990), and Moss (1992). De- 
tailed calculations of the phase mixing process have been 
made by Charboneau & McGregor (1993) and Sakurai et 
al. (1995). 

3.3. Stronger initial fields, and interim conclusions 

If Bq > B\ (Eq. |2^), torques become effective before rota- 
tional smoothing can make the field axisymmetric. If the 
initial field is non-axisymmetric, it will therefore stay non- 
axisymmetric. Such a field still has magnetic surfaces how- 
ever. Alfven waves and the slow mode continuum travel 
on these surfaces (Goossens et al. 1985), causing phase 
mixing just as in the axisymmetric case. These waves will 
damp out on the time scale t p . The final field will then be 
similar to the initial field. But what is the state of rotation 
at the end of the process? Since the magnetic surfaces are 
not axisymmetric, and differential motions are damped on 
these surfaces, the stationary state is one of uniform rota- 
tion, and nonuniform rotation is possible only if the initial 
field is axisymmetric (Mestel et al. 1988, 1990). 

This completes the discussion for the case when only 
the azimuthal equation of motion is taken into account. 
Depending on the initial field strength, the final state in 
this picture is either a uniformly rotating magnetic star 
with a non-axisymmetric field, or a differentially rotating 
star with an axisymmetric poloidal field, with rotation 
constant on magnetic surfaces. This is only a preliminary 
scenario, of course, since we have allowed only for purely 
azimuthal motions. The axisymmetric poloidal end-state, 
if it ever were to materialize, would evolve further by mag- 
netic instabilities. 
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4. Magnetic instabilities 

The fluid motions, taken to be purely azimuthal in the 
above, are now unconstrained, and the evolution of the 
field becomes much more interesting since magnetic insta- 
bilities can take place. 

Consider first the case where the initial field is weak, 
such that differential rotation has time to produce a pre- 
dominantly azimuthal field by rotational smoothing (Sect. 
3.1). Stronger initial fields are taken up again in Sect. 0. 
The restriction to weak initial fields isolates the magnetic 
instabilities that appear from the phase mixing process 
discussed above, which can be ignored here. 



4-1- Magnetic shear instability 

In a differentially rotating object with angular velocity 
decreasing outward, a weak field is generically unstable 
to magnetorotational instability (Velikhov 1959, Chan- 
drasekhar 1960, Fricke 1969, Acheson 1978, Balbus & 
Hawley 1991, 1992, Balbus 1995). This is a linear insta- 
bility in which the magnetic field mediates the release of 
free energy in differential rotation. The growing magnetic 
perturbations couple the fluid at different radii, the dif- 
ferential rotation acting on them converts the free energy 
into magnetic and kinetic energy on smaller scales. The 
instability acts when the rotation rate decreases, but an- 
gular momentum increases outward. In other cases, i.e. 
when angular momentum (</) decreases, or angular veloc- 
ity increases outward, ordinary hydrodynamic shear insta- 
bilities are known to exist. These cases are less relevant 
for differentially rotating stars or accretion disks, however. 
Whether purely hydrodynamic instabilities exist for f2 de- 
creasing and J increasing with radius is still controversial. 
The possibility that this case might actually be hydrody- 
namically stable should be taken seriously, however, given 
that instabilities have so far been found neither theoreti- 
cally, nor numerically, nor experimentally. A physical ar- 
gument for stability, not constituting a proof, has been 
given by Balbus et al. (1996). 

Where magnetic shear instability is present, it gener- 
ates a magnetic form of turbulence. In stars, a strongly 
stabilizing factor is the stratification. Except in regions 
close to convective instability, this limits the instability to 
cases with strong differential rotation, and a rotation rate 
near the maximum value. Apart from these extreme sit- 
uations, magnetic shear instability is limited to displace- 
ments on horizontal surfaces, and redistributes angular 
momentum over such a surface if the rotation rate de- 
creases with cylindrical radius (Kato 1992, Balbus 1995). 
The effect of the instability is thus much more limited 
than in the case of an accretion disk, where differential 
rotation is dominant, pressure effects small, and magnetic 
shear instability endemic. 

On sufficiently small length scales, however, thermal 
diffusion can strongly reduce the stabilizing buoyancy ef- 



fect of the thermal stratification. As in the case of purely 
hydrodynamic instabilities (Zahn 1974, 1983), one would 
therefore expect that i) linear magnetic shear instability 
should reappear at large wavenumbers when thermal dif- 
fusion is taken into account, and ii) angular momentum 
transport by some form of small scale magnetic turbu- 
lence will take place. The first of these was already demon- 
strated to be the case in the classical and detailed analysis 
by Acheson (1978). He showed that the necessary and suf- 
ficient condition for linear instability is 

N 2 n 

g=-dlnr!/dlnr> — (29) 

(cf. eq [[3]), where N is the buoyancy frequency and k 
the thermal diffusivity. For conditions in stellar interiors, 
K 3> n. This is because the magnetic diffusion involves the 
random walk of electrons, which have much shorter mean 
free paths than the photons that carry the heat. Thus 
magnetic shear instability is possible with weak rotation 
gradients |dlnf2/dlnr| <C 1, provided that the rotation 
rate is not too slow compared with the buoyancy fre- 
quency. For a star rotating as slowly as the Sun, however, a 
fairly strong gradient in rotation rate is still needed for in- 
stability. For the present Sun, with N ~ 1CT 3 , = 3 1(T 6 , 
tj/k ~ 5 10 -5 , instability would require q ~ 0(1). This is 
a much larger degree of differential rotation than allowed 
by current observations (Schou et al. 1998) of the Sun's 
interior rotation. 

How much angular momentum transport this process 
is likely to produce, in cases where a large enough shear 
exist, is not certain at the moment. The issue is discussed 
in Spruit & Balbus (in prep.). Several lines of argument 
presented there suggest that a modestly effective trans- 
port of angular momentum is possible. A simple, perhaps 
convincing, argument starts by assuming that the mag- 
netic turbulence has the effect of enhancing the magnetic 
diffusion of large scale fields. The amplitude at which this 
turbulence saturates is then found by requiring that insta- 
bility is just possible according to condition (p9|), with the 
magnetic diffusivity w replaced by an effective diffusivity 
n c . This condition then determines n c : 



rj e = 2qn 2 K /N 2 . 



(30) 



[Recall that the same argument, applied to convective in- 
stability, can be used to derive the familiar mixing-length 
model of convection.] The fluid motions in the turbulence 
are of slow-mode and Alfvenic type, so that the kinetic 
energy is of the same order as the energy in the magnetic 
field. If v c is the effective viscosity due to the fluid motions, 
the magnetic Prandtl number P c = v c jr] c is then of order 
unity. Numerical simulations suggest P c « 0.1 (Hawley et 
al. 1995, Brandenburg et al. 1995, Matsumoto & Tajima 
1995). With the assumptions made, ( ^9| ) now yields the 
effective viscosity expected: 



n 2 

v c = 2P c q K —. 



(31) 
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Since the induction equation is linear, the winding- 
up, rotational smoothing and phase mixing processes dis- 
cussed above all just proceed as before, the only differ- 
ence being that the magnetic diffusion has to be replaced 
by the effective value (|3(]). Quantitatively, the difference 
is not very large. For the initial Sun assumed above, the 
effective diffusivity in the presence of magnetic shear in- 
stability would have been only 10-100 times larger than 
the microscopic value. With this effective diffusivity, the 
rotational smoothing time (Q) decreases only by a factor 
of order unity. 

The net effect of magnetic shear instability therefore 
may be small. Even though it is not always suppressed by 
a stable thermal stratification, its effects are modest. This 
conclusion, however, is preliminary, as it depends on an 
estimate of its nonlinear development. Numerical simula- 
tions should be able to settle this question more defini- 
tively. 

4-2. Instabilities of an azimuthal magnetic field 

During the rotational smoothing and phase mixing pro- 
cesses in an initially weak field, the azimuthal field is much 
larger than the poloidal component. An azimuthal field 
has its own instabilities. Magnetic shear instability oper- 
ates also on such an azimuthal field, but is more properly 
regarded as a form of shear instability, extracting its en- 
ergy from differential rotation (as discussed above) . Other 
forms of instability exist, which draw their energy from the 
magnetic field itself. They are of two distinct types. Per- 
haps the most well-known is magnetic buoyancy or Parker 
instability (Parker 1966), due to a vertic al gr adient in the 
field strength. This is discussed in Sect. |5.3| . The second 
is a pinch-type instability. 

5. Tayler instability 

The second form of instability derives its energy from 
(nearly) horizontal interchanges. The classical result in 
this context is that of Tayler (1973), who considered adi- 
abatic perturbations (rj = n = 0) in a stratified, nonro- 
tating star. By an energy method he showed that every 
purely azimuthal field B^{r, 9) in a stably stratified (non- 
rotating) star is unstable on an Alfven time scale, no mat- 
ter how weak the magnetic field. This form of instability 
is closely related to many forms of instability in pinches 
and torus configurations (e.g. Tayler 1957). Its behavior 
under the strongly stratified conditions in stellar interiors 
is sufficiently distinct, however, that I find it convenient to 
refer to it by the separate name of 'Tayler instability'. In 
poloidal field configurations, pinch-type instabilities also 
occur in stellar interiors; they were studied by the same 
method by Wright (1973), see also Sect. [| 



5.1. Previous results 

The most unstable motions are to = and m = 1 modes 
with nearly horizontal displacements. The displacements 
take place at nearly constant total pressure, i.e. the Eule- 
rian perturbation S(P + B 2 /8tt) vanishes. They are essen- 
tially local in the r and directions: when the instability 
condition is satisfied at a point (r, 6) , unstable perturba- 
tions can be found that are confined to a small neighbor- 
hood of this point. In spherical coordinates (Goossens et 
al. 1981; Tayler's was derived in cylindrical coordinates), 
the necessary and sufficient conditions for instability are 

d e \n{B 2 sin 6 cos 9) > 0, (to = 1) (32) 

cos 9 d e MB 2 / sin 2 9) > 0, (m = 0) (33) 

From the induction equation it is evident that azimuthal 
fields created by winding up of a poloidal field must have 

B ~ 6 (6 < 1) (34) 

near the pole. Hence there is always a region near the pole 
where condition ( |32"| ) is satisfied, and the field is unstable. 
Away from the pole, the most unstable modes can be ei- 
ther to = 1 or to = 0. As an example, consider the field 
resulting from the winding up of an initially uniform field 
parallel to the rotation axis, when rotation is a function 
of r only. This field has 

B^ ~ sin 9 cos 9. (35) 

The region 7r/4 < 9 < 37r/4 is stable for to = 1 modes, 
the polar caps are unstable, and the axisymmetric modes 
are stable. The growth rate is of the order (Tayler 1973, 
Goossens et al. 1981) 

cr w n A = v A /r (ft < n A ). (36) 

The unstable displacements are sketched in Fig. [|. 

When rotation is included, the energy method does not 
work, and somewhat weaker results are obtained. Never- 
theless, Pitts & Tayler (1986) showed that rotation does 
not by itself remove the instability. There are still unsta- 
ble to = 1 modes, though their growth rate is reduced. 
When rotation is rapid compared with the Alfven travel 
time, as I have assumed above, the growth rate found by 
Pitts and Tayler is of the order 

a k. Q 2 A /Q (fl > n A ). (37) 

This dependence on the rotation rate is typical for in- 
stabilities of an interchange type. It represents the fact 
that the Coriolis force, being perpendicular to the veloc- 
ity, does not enter into the energy budget of the perturba- 
tions, though it can still affect the shape of the unstable 
modes and their growth rates. 

The growth rate ( |37j ) is still fast compared with the 
rotational smoothing and phase mixing time scales, so 
that the effects of the instability are likely to significantly 
change the evolution of an initially weak field as sketched 
in Sect. 0. 
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Fig. 1. Unstable displacements in an azimuthal field near the 
pole. Shown is the m = 1 mode, which occurs under the widest 
range of conditions. The displacements are along horizontal 
surfaces (indicated by arcs) 



5. 2. Effects of rotation and diffusion 

The published results deal mostly with the stability and 
growth rates considering the star as a whole, using an 
energy method. Tayler's result shows, however, that the 
essence of the instability is a local interchange-type pro- 
cess (see Appendix 2 in Tayler 1973). Hence it is natural 
to study the instability by a modal analysis, in a local 
approximation. This also has the advantage that the ef- 
fects of rotation, magnetic and thermal diffusion can be 
included, which is not possible with an energy method. 

The result of such a local approximation, for a purely 
azimuthal field, has been given by Acheson (1978). His dis- 
persion relation includes the effects of rotation, viscosity, 
thermal and magnetic diffusion, and is rather complicated. 
The special case of an azimuthal field at the equator has 
been analysed in detail by Acheson (1978, 1979). Tayler 
instability, however, disappears at the equator, and shows 
its most characteristic behavior at the poles. 

In the following I derive the properties of the insta- 
bility by heuristic arguments. In the Appendix, I derive 
these results more systematically from the dispersion rela- 
tion. I also show there that Acheson's local approximation 
agrees, in cases where a comparison is possible, with all 
results based on the more rigorous energy method. This 
establishes the validity of his local analysis for the case of 
Tayler instability. The Appendix also discusses a few sub- 
tleties that are not captured by the heuristic derivations. 

A useful simplification is to ignore viscosity. In a stel- 
lar interior in which radiation is not the dominant source 
of viscosity, it is smaller than the next larger diffusion 
process, magnetic diffusion, by 1-2 orders of magnitudes. 
Use cylindrical coordinates (ro, <p, z) at the pole, where the 



1 The effects of viscosity may need to be reconsidered, how- 
ever, for the interiors of massive stars, where radiation pressure 
is important. 



rotation axis is in the z direction. The magnetic field is 
B = (0, B(zu), 0). Define an Alfven frequency wa: 



B 



(4 7 rp) 1 /2 T 



(38) 



(which is a function of zu, as opposed to the very similar 
quantity f^A used above). Since the analysis is local, the 
field is characterized completely by two numbers, the local 
value B of the field strength, and its radial gradient 



p — d In B/d him. 



(39) 



The vertical gradient of the field is neglected here. Insta- 
bilities associated with such a gradient are discussed below 
in Sect. 



5.3 



The rotation rate is assumed to be uniform 
in the present analysis. A gradient dQ/dz is present and 
causes its own forms of instability, but these have been 
covered already above in the discussion of magnetic shear 
instability. 

There are several different time scales involved in the 
problem. The largest frequency is the buoyancy frequency 
N, reflecting the strong effect of stratification. The next 
lower frequency is the rotation rate f2. For the weak fields 
that turn out to be most relevant, the Alfven frequency wa 
is small compared to f2. In this case we have the ordering 



n > n > uj a . 

The perturbations are of the form 



(40) 



(41) 



and are essentially incompressiveel, divv = 0. In the ab- 
sence of diffusion (re = r\ = 0), the growth rate of the 
instability is maximized for n — > (Pitts & Tayler 1986). 
In this limit, the instability is essentially confined to hor- 
izontal surfaces. The growth rate, which in the absence 
of rotation is of the order wa, is reduced by the Coriolis 
forces to a value of the order 



(42) 



but the instability condition itself does not depend on f2. 
This behavior is generic for interchange-type instabilities 
in the presence of rotation, as discussed above. The insta- 
bility conditions are (Tayler 1957, see Appendix) 



P> 



1 (rn^O), 



and 



p > 1 (to = 0). (43) 



The m = 1 mode thus occurs under the widest range 
of conditions (usually, its growth rate is also largest if 
several modes are unstable). In the following, I restrict 
the discussion to to = 1. 



2 Small expansions and contractions occur in order to make 
the total pressure perturbation vanish for the most unstable 
modes. These disappear in the limit va/c s ~ lua/N <C 1. 
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5.2.1. Magnetic diffusion only 

Consider first the case where magnetic diffusion is in- 
cluded but thermal diffusion is neglected. This would be 
appropriate for situations where the dominant contribu- 
tion to the buoyancy frequency is a composition gradient. 
In order to avoid doing work against the stable stratifica- 
tion, the unstable displacements must be nearly horizon- 
tal, v z /v p »j l/n <C 1. For displacements of amplitude £, 
the work done per unit mass against the stable stratifi- 
cation is \^ 2 {l/n) 2 N 2 . The energy gained from the field 
configuration is \u 2 pj; 2 . For instability, the field must be 
strong enough, such that oj\ > (I 2 /n 2 )N 2 . For a given 
vertical wavenumber, this is most easily satisfied at the 
longest possible horizontal wavelength, which is of the 
order r, the (spherical) radius. The vertical wavenumber 
thus has to satisfy 



n 2 > 



N 2 



where tt = (?t. 2 k) _1 is the thermal diffusion time at 
wavenumber n, and r the time scale of the process un- 
der consideration. In our case r is the adiabatic instability 
time scale c -1 (Eq. [12]). For these time scales, the temper- 
ature perturbations are reduced by the factor in brackets 
in (pf). Since k 3> rj, this factor is approximately t/tt- 
The same argument as that leading to ( f47j ) then yields as 
conditions for instability (cf. Appendix eq |A2E ) p > — h 
and 



UJA 

n 



> 



(44) 



With this condition the critical field for instability in 
the present Sun is now about 100G, and the typical growth 
time (from|4^) is of the order 10 5 yr (for conditions signifi- 
cantly exceeding marginal) . The wavenumber at marginal 
stability is nr = 10 3 , corresponding to a wavelength of 
3000km. 



At large wavenumbers, however, diffusion starts affecting 
the perturbations. The condition that the rate at which 
they decay by magnetic diffusion does not exceed the 
growth rate yields 



n i] < a, 



(45) 



where a is the growth rate in the absence of stratification 
and diffusion. For our ordering Q 3> lu a , a is given by 
(fl2"|). Wavenumbers for which both conditions ( j44| , f45| ) are 
satisfied then exist if 2 conditions are satisfied. The first 
is 



1 



(46) 



which is always satisfied in some region near the poles. 
The second is 



LO A > fi(^) 1 /2(JL)l/4. (n » K = 0) 



(47) 



This condition applies to the most unstable mode, m = 1 
. For example, in the outer core of the present Sun (r ~ 
5 10 10 , N ~ 10~ 3 , 7] ~ 2 10 3 , k ~ 4 10 7 , ft = 3 10~ 6 ), the 
minimum azimuthal field strength needed for instability, 
B = uj A r(Airp) 1 / 2 , would be of the order 1000G. 

5.2.2. Magnetic and thermal diffusion 

Condition ( |47| ) overestimates the field strength required 
for instability if the stratification is due to a thermal gra- 
dient. At the short wavelengths where magnetic diffusion 
plays a role, thermal diffusion strongly reduces the sta- 
bilizing temperature perturbations. As in the case of hy- 
drodynamic instabilities in a stable thermal stratification 
(Zahn 1974), this effect can be taken into account by re- 
placing N with an effective buoyancy frequency N, 



N 2 = N 2 /(1 + t/t t ), 



(48) 



5.2.3. Effect of a composition gradient 

The situation is more complicated if the stratification is 
due to the combined effects of a thermal gradient and a 
gradient in composition. If \x is the mean atomic weight 
per particle, the buoyancy frequency is given by (e.g. Kip- 
penhahn & Weigert, 1990): 



N 2 =N^ + N 2 , 



(50) 



where the thermal and compositional contributions 
and are given by 



#r = ^(V-V a ); N 2 = ^V fl 



(51) 



with 

„ dlnT „ 
dlnP' 



dlnT 
d\nP 



dhip 
dlnT 



p 4 , 



91np s 
d In n j p rp 



(53) 



Here straight derivatives measure the variation of the 
physical variables with depth in the stratification, partials 
are thermodynamic derivatives of the equation of state, 
and S is the entropy. For an ideal gas, (P = plZTf/j,), 
<j> = 6 = 1, hence 



iV| = |(V 



V a ); 



(54) 



For adiabatic perturbations, the buoyancy frequency in 
the formulas above is just replaced by ((5^). For nonadi- 
abatic perturbations, we need to take into account that 
inhomogeneities in composition and temperature diffuse 
at different rates. The diffusivity of inhomogeneities in 
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ji is of the same order as the viscosity, since both quanti- 
ties scale with the mean free path and velocity of the ions. 
We can then neglect compared with the magnetic dif- 
fusivity. The same arguments as those leading to eq ( [i"7j ) 
and ( |49| ) then gives the instability conditions p > — | and 



Hence we can concentrate on here the case H <C r. 
In terms of the Alfven frequency, the instability condition 
can be written as 



7 V H , r9 
> — -—N 2 . 

—p k r 



(58) 



UJA 

n 



> 




(55) 



(ft > OJA, («>!)) 

5.3. Magnetic buoyancy instability 

If the magnetic field strength increases in the direction of 
gravity, the gas is supported against gravity in part by the 
magnetic pressure. If the gradient is sufficiently strong, the 
free energy in the field gradient can be released by buckling 
of the field lines. As with all instabilities where displace- 
ments in the vertical are necessary, the stratification pro- 
vides a strong stabilizing force against such buckling. Like 
in the other instabilities, however, this stabilizing force is 
reduced by thermal diffusion on sufficiently small scales. 
To demonstrate the properties of the instability it is suf- 
ficient to consider the situation at the equator. This case 
has been analysed in detail by Acheson (1978). Ignoring 
viscosity, the condition for instability is (Acheson's Eq. 
7.27): 



ft 2 din ft 



din 7' 



-'if- 2 



dlnB n-/N 2 



dlnr 



K W 



(56) 



where uja — Va/t, H the pressure scale height, and 7 the 
ratio of specific heats. The first term d escri bes magnetic 
shear instability, discussed above in Sect. 4T. By assuming 
uniform rotation, we get: 



r 
H 



2)p> 



rj -fN 2 



(57) 



where p = dlni?/dlnr. This condition is independent of 
the rotation rate. The growth rates, however, are reduced 
by rapid rotation. Compared with the nonrotating case, 
they are smaller by a factor cr/ft, where a is the growth 
rate in the absence of rotation. 

Near the center of the star, r — ► and H — > 00, so that 
the second term in the bracket dominates. Instability is 
then possible if the field strength increases outward. This 
is obviously not the normal buoyancy instability. It is, in 
fact, the same pinch- type instability as the instability that 
appears near the poles, and is associated there with the 
horizontal gradient of the field strength. These have been 
discussed in Sect. |5[ The directions of the rotation axis 
and stratification are different in the present case, but in 
the absence of rotation and gravity the instability would 
be the same. 



For a smooth field gradient, —p ~ 0(1), the critical Alfven 
frequency for instability is of the order (subscript b for 
buoyancy) : 



^Ab 



C^y/ 2 N. 



(59) 



Comparing this with condition (|4£j) for diffusive Tayler 
instability, with its critical Alfven frequency a; at, we have 



^Ab _ ^Y /2 {-) 1/A { — ) 1/2 (^—"> 1,A 
ft k r 



LUAT 



I] 



(60) 



This is minimized for the largest possible rotation rate, 
ft ~ N. The dominant factor is the last one, which is quite 
large since the magnetic diffusivity is so small for stellar 
length and time scales. One finds that u>Ab / ^at 3> 1 for 
main sequence stars, white dwarfs and neutron stars. For 
the present Sun, uj Ab / ^at ~ 10 3 . 

Buoyancy instability thus requires much stronger fields 
than Tayler instability. A field wound up by differential ro- 
tation into an azimuthal field therefore becomes unstable 
to Tayler instability first. If the instability is able to limit 
the growth of the toroidal field, a slowly wound-up field 
will settle at a value near the marginal conditions for in- 
stability (Mestel & Weiss 1987), and it is unlikely that 
buoyancy instability will become important. 

The stability conditions discussed are summarized in 
Fig-! 

6. Poloidal field instabilities 

Though azimuthal fields are conceptually attractive as a 
natural result of differential rotation, the possibility of 
purely poloidal fields (B^ = 0) or general mixed poloidal- 
toroidal fields must also be considered as possible initial 
field configurations. As discussed further in Sect. 0, few re- 
sults exist on the stability of mixed poloidal-toroidal fields. 

For purely poloidal fields, strong results are again 
available. The most important result is that, most likely, 
all purely poloidal fields in stars are unstable to adiabatic 
perturbations, in the absence of rotation. This was demon- 
strated by Wright (1973) and by Markey & Tayler (1973, 
1974). These authors considered poloidal fields in which 
some or all field lines in each meridional plane are closed 
within the star. Then on each of these planes there is (at 
least one) point where the field strength vanishes. Near 
this point the field lines are ellipses centered around the 
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Fig. 2. Stability conditions for an azimuthal magnetic field 
as functions of field strength and rotation rate, uja is the 
Alfven frequency Va/t, N the buoyancy frequency, q is the dif- 
ferential rotation d In fid In r. Assumed stellar parameters are 
77/jt = 10 -4 , r 2 N/ri = 10 15 . With increasing field strength, the 
first instability to appear is Tayler instability. Magnetic shear 
instability appears above a minimum rotation rate. 

point. The configuration thus closely resembles the con- 
figuration of an azimuthal field near the pole, and one ex- 
pects similar pinch-time instabilities, driven by the curva- 
ture of the field linescl. The situation differs in the present 
case in that the direction of the stable stratification is 
now within the plane of the field lines. Unstable displace- 
ments £ must be nearly incompressive and close to hori- 
zontal surfaces, g-£ = 0. Displacements in latitude, with a 
small length scale in the azimuthal direction satisfy these 
requirements. They have the same effect as an m = 1 
displacement in an azimuthal field near the pole. Wright 
(1973) and Markey & Tayer (1973, 1974) find that these 
are indeed the most unstable ones, and that they make all 
poloidal fields with closed field lines unstable. The growth 
rate, as expected, is of the order lo a - 

A special case occurs when none of the field lines is 
closed inside the star. All field lines then cross the stel- 
lar surface. An example would be a uniform field inside 
the star, with a dipolar vacuum field outside. The results 
by Markey and Tayler do not apply to this case, but a 
simple argument (Flowers & Ruderman 1977) shows that 
this case is equally unstable. As a trial function for the dis- 
placements consider splitting the star in half by a plane 
containing the axis, and rotating one of the halves over 
180°. This changes neither the thermal, nor the gravita- 
tional, nor the magnetic energy of the star. It changes the 
external vacuum field, however. By a suitable choice of the 
plane, the rotation brings fields of opposite polarity closer 



Note that this is now the third instance of such instabilities. 
The first time was instability at the poles, the second time in 
the instability of an azimuthal field near the center of a star in 
Sect. |J. 



together on the surface, which lowers the energy of the 
external vacuum field. This is most easily visualized for a 
dipole field. The star is then analogous to a bundle of bar 
magnets, oriented in parallel. By splitting the bundle in 
two and rotating one half, the energy of the bar magnets 
in the field of the other half is reduced by an amount of 
the order B 2 V/8ir, where V is the volume of the star. Note 
that this energy change is due entirely to the external field 
energy, the internal field energy does not change. With this 
value for the change in energy, the growth rate of the insta- 
bility is of the order wa- The process can be continued by 
rotations along other planes, which increases the complex- 
ity of the external field (higher multipoles), and reduces 
its strength. Asymptotically, the external field energy and 
the growth rate of the instability vanish. 

6.1. Effect of rotation 

For adiabatic perturbations, the effect of rotation has been 
studied by Pitts & Tayler (1986). The results are less 
complete than in the nonrotating case, since the power- 
ful energy method fails for rotating systems. For the cases 
studied, however, poloidal fields were again found to be 
unstable. The reason for this can be visualized easily by 
considering the most unstable displacements of the non- 
rotating case. These have azimuthal wavenumber m 3> 1 
and £e ^ £,r,(,<f>- The Coriolis force on such displacements 
is in the azimuthal direction. On neighboring meridional 
planes, such that the phase of the perturbation differs by 
7T, the Coriolis forces are opposite. These forces can be bal- 
anced entirely by the azimuthal pressure perturbation, so 
that the net effect of the Coriolis force vanishes for high- 
m perturbations. The instability then proceeds under the 
same conditions and with the same growth rate as in the 
nonrotating case. An explicit example of this effect has 
been given in Spruit & Taam(1990), where poloidal field 
instability was studied in the very analogous case of a uni- 
formly rotating disk. Numerical simulations of this insta- 
bility in rotating disks have been made by Stehle (1998). 

One should therefore expect rotation to have even less 
effect on poloidal fields than it has on the azimuthal fields 
discussed in Sect. ||. 

6.2. Effect of magnetic and thermal diffusion 

The effect of the diffusivities has not been studied for 
poloidal fields. Given the nature of the unstable diplace- 
ments, which are closely analogous to the m = 1 instabil- 
ities near the pole in an azimuthal field, we should expect 
the effects of the diffusivities to be qualitatively the same. 
In the nonrotating azimuthal field case, diffusion does not 
stabilize the configuration at any field strength, but only 
affects the growth rates. Since the effect of rotation on 
the unstable displacements is small for poloidal fields (see 
above), it is a good guess that poloidal fields will be unsta- 
ble at all field strengths, even in the presence of rotation 
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and magnetic diffusion. Demonstration of this may need 
more detailed study, however. 

7. Initial field strong: stable equilibria? 

In the above, I have assumed the initial field to be weak, 
so that a predominantly azimuthal field quickly develops 
by differential rotation. If the initial field is strong (as 
measured by condition ^2j), this is not the case. Instead, 
in this case, the differential rotation can be regarded as 
a perturbation of the field which damps out on the phase 
mixing time scale (p7j). In doing so, the field settles to a 
stable equilibrium state, if one exists. 

This is a loose end in the story since in spite of exten- 
sive work (mainly in the 60's and 70's, see e.g. Tayler 1980, 
Borra et al. 1982, Mestel 1984), it is still not known how 
to prove the existence or absence of stable magnetic equi- 
libria in stars, whether they rotate or not. As discussed 
above, all purely toroidal field configurations are unsta- 
ble above a critical field strength given by the diffusive 
Tayler instability condition ( |49| ) . All purely poloidal fields 
are likely to be unstable as well, as discussed in Sect. [| 
For mixed poloidal-toroidal fields, stability analyses ex- 
ist as well (e.g. Wright, 1973, Tayler 1980), but they do 
not lead to very general conclusions. As argued by Mestel 
(1984) and Tayler (1980) this leaves the possibility that 
special stable configurations might exist with poloidal and 
toroidal fields of similar strength, but no example of such 
a field is known. 

An indirect, and not completely compelling, argument 
that stable configurations are possible can be made by ap- 
pealing to the conservation of magnetic helicity (Moffat, 
1989, private communication). If A is the vector potential 
of B (with a suitable gauge) , then it can be shown that the 
integral of H = B • A over volume is a conserved quantity 
in ideal MHD (i.e. in the absence of magnetic diffusion) 
(Woltjer 1958, see also Taylor, 1974). If a star is born with 
a field for which H is nonzero, then this field must reach 
a stable equilibrium configuration. Suppose it is not ini- 
tially in a stable equilibrium. By putting in a sufficient 
damping mechanism such as a viscosity, we can make sure 
that the energy released by the instability is dissipated on 
some finite time scale. The final state must then be a stable 
equilibrium at a finite field strength, since an infinitesimal 
field can not have a finite helicity. This argument is unfor- 
tunately not compelling, since it is possible that unstable 
magnetic fields can evolve in such a way as to develop sin- 
gularities (current sheets) within an effectively finite time. 
Once such current sheets develop, reconnection sets in and 
the helicity is no longer conserved. 

While theory fails to give a clear answer to this long 
standing question (Mestel 1984), we may appeal to obser- 
vations to argue that stable configurations do in fact exist 
in stars. The magnetic A stars and the magnetic white 
dwarfs have strong magnetic fields (of the order 10 4 and 
10 7 G, respectively), that do not change on time scales of 



at least decades. Since these stars are also slow rotators, 
the field would change on an Alfven time scale, of the order 
of a year and a day, respectively, if their field configura- 
tions were not stable. The A stars have convective cores, 
in which a steady dynamo might possibly exist to pro- 
duce the observed stable field. The magnetic white dwarfs 
do not have a plausible location for such a dynamo. Very 
cool white dwarfs may form a crystalline lattice in part 
of the interior, which might be able to anchor a strong 
field. The known magnetic white dwarfs (Schmidt & Nor- 
sworthy 1991, Liebert 1995) however, are not all of this 
type. We can therefore take the magnetic white dwarfs as 
a fairly strong argument for the existence of stable field 
configurations in stars. 

Assuming that this is the case, initial field configu- 
rations with strength above (^2|) will evolve into stable 
configurations on the phase mixing time scale. The final 
field strength would depend on the degree to which re- 
connection has taken place during the evolution of the 
configuration. The star would be uniformly rotating, ex- 
cept for special axisymmetric field configurations aligned 
with the rotation axis (which are apparently not realized 
in the observed magnetic A stars and white dwarfs). 

8. Discussion 

In the above I have reviewed the known processes relevant 
for magnetic fields in differentially rotating stably strati- 
fied stellar interiors. This excludes dynamo processes such 
as are thought to occur in convective zones. In initially 
weak fields, (condition |l|, about 30G for the Sun) the 
non-axisymmetric components are smoothed out by ro- 
tational expulsion, producing an axisymmetric field. Dif- 
ferential rotation is subject to damping by phase mixing. 
This results in a uniformly rotating star if sufficient time 
is available, compared with time scale on which the inter- 
nal structure evolves and compared with the spin down 
timescale due to external torques. 

The azimuthal field that results from winding-up of an 
initially weak field is subject to instabilities. In the stably 
stratified environment of a stellar interior they are of two 
types. There are the Parker (or magnetic buoyancy) and 
Tayler (or stratification-modified pinch-type) instabilities, 
both driven by the magnetic field energy in the toroidal 
field. They are of particular relevance for fields produced 
by differential rotation, through the winding-up of a weak 
initial field. 

In addition to these instabilities driven by the free 
energy in the magnetic field, there is magnetic shear in- 
stability, which feeds on differential rotation. It occurs in 
any field configuration, and occurs already in very weak 
fields. Due to effects of stratification and magnetic diffu- 
sion, however, the rotation gradient has to be significant 
for it to occur, especially in slowy rotating stars. Poloidal 
fields configurations have similar instabilities. 
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It turns out that Tayler instability (pinch-type insta- 
bility in the presence of a stabilizing stratification) is of 
particular relevance. In an azimuthal field, it generically 
occurs in a region near the pole, in the form of an m — 1 
displacement of the field lines along horizontal surfaces. 
This instability is of a local interchange type, and the 
effects of rapid rotation and magnetic and thermal diffu- 
sion can be included in the analysis. I find that Acheson's 
(1978) dispersion relation can be applied to this form of 
instability. The results of Sect. || and the Appendix show 
that Tayler instabilities probably are more relevant than 
the better known buoyancy (or Parker-) instabilities. In 
stars with magnetic fields that have been wound up by dif- 
ferential rotation, they set in at the lowest field strength. 

One can wonder how complete our 'catalog' of known 
instabilities is. As discussed above, the situation for nearly 
azimuthal fields, such as would result from differential ro- 
tation, is probably quite satisfactory. The energy method 
used by Tayler (1973) gives necessary and sufficient con- 
ditions for the adiabatic, nonrotating case, and is thus 
complete for this case. Quite important is that this anal- 
ysis shows that the instabilities are of a local nature (in 
the r and 9 coordinates). This allows a complete study 
of the effects of rotation, viscosity magnetic and thermal 
diffusion by Acheson's (1978) approach. Barring the possi- 
bility of diffusive non-local instabilities that don't have an 
equivalent in the nonrotating adiabiatic case, the stability 
of azimuthal fields can thus be analyzed completely. The 
same applies essentially to purely poloidal fields, where the 
energy method again shows the instabilities to be local. A 
complete study of the effects of rotation and diffusion on 
these instabilities still has to be done for this case, how- 
ever. 

Another question is, of course, how the magnetic con- 
figuration would evolve nonlinearly under these instabil- 
ities, in particular how effective they would be at trans- 
porting angular momentum. It is conceivable that signifi- 
cant progress in this question can be made with numerical 
MHD simulations. 

The stability of initially strong fields, such that the 
differential rotation is not strong enough to wind the field 
into an azimuthal configuration, is as open a question as 
before. The same applies to cases where the field has gone 
through phases of winding-up and phase mixing. In the 
absence of forces that continue to create differential rota- 
tion, the phase mixing process eliminates differential ro- 
tation, and leaves the field in a configuration of unknown 
stability. If such a field is to be stable, it is clear that it 
can neither be a purely poloidal nor a purely toroidal field 
since these are unstable on short time scales. The possi- 
bility that a mixed poloidal-toroidal configuration can be 
stable on long time scales can not be excluded (Tayler 
1980, Mestel 1984). The magnetic white dwarfs are fairly 
strong observational evidence that such configurations do 
in fact exist. 
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A. 



Appendix: 
diffusion 



Tayler instability with rotation and 



n 2 + I 2 



lu a = mV A /-cu, 



dh =Qvj~ (l/n)d z , 



lP/p, 



P 

E = In—, 

P 7 



F = ln— , 

pw 



P = 



dlnB 
dlnnj' 



(A3) 



(A4) 



(A5) 



and H the pressure scale height. The magnetic field 
strengths relevant here are weak in the sense Fa € c s , 
i.e. the magnetic pressure is negligible compared to the gas 
pressure. The terms of order v A /c s on the third and fourth 
lines can then be neglected, and the gradient of the den- 
sity in the w-direction can be neglected. Under the same 
assumption, the expected growth rates, of the order lua or 
smaller, are small compared with the buoyancy frequency 
N, hence instability is possible only for l/n ~ lua/N <C 1. 
The terms with factor I jn in the first and fourth lines can 
thus be neglected. The same assumptions also imply that 
s/n ph 1. Thus we have 



vjdh.F ~ dln(S/ro)dlni<7 = p — 1, 



(A6) 



We start with the local dispersion relation in Acheson 
(1978). The goals of the analysis are to show that this 
relation yields the correct results for Tayler instability in 
the known cases, and to extend these to include rotation, 
magnetic and thermal diffusion. 

The perturbations are of the form 



a i(lm-\-mtp-\-nz~ ut) 



(Al) 



The field is axisymmetric and purely azimuthal. Its 
strength B is taken to be a function of cylindrical radius 
w only. If B depends on z as well, additional forms of 
instability are possible (buoyancy or Parker instability). 
This is discussed in Sect. 5.3. At the pole of a uniformly 



rotating star, where gravity is parallel to the rotation axis, 
and ignoring viscosity, the dispersion relation reduces to 



Yl 

w 



+ 



2fim 



I W LU 

2 + — 



IKS 



n H <x>7 + ins 1 



d h F 



to 



d h E^—d z HP/ P ~<) 

n 



I iV 2 7 
n g 



(A7) 



The term involving E on the first line can be neglected, 
but not the term on the second line, since it is multiplied 
by g. For u ~ lua, this term is of the same order as the 



first term on the second line. Eq. (A2) thus reduces to 



V 2 

-4(2nm + 2w)(p- 1) 
w 



N 2 



2n 



u + ir\s n- 



1 U! + IKS 2 /j 



\u)(u + irjs 2 ) — uPp 



1) 



2n(uj + irjs 2 ) + 2 



ml^ 1 



0.(A8) 



W7 + IKS' 



irjs' 



I 



n W7 + Iks 



vl 



Cg LU + IKS 2 /7 



2n 



, A . a (p- 



\u)(u> + irjs ) 



1) 



;dhE To see that the stability conditions implied by this rela- 
tion are the same, in the adiabatic case as, those found by 
Tayler (1973), Goossens et al. (1981), and Pitts & Tayler 
1(1986), set 77 = k = 0, so that 



VI 



I 2 . 
~ N 



(- 2 ^i) 



2il 



, V? U) + IKS 2 

ir]S + : — tti-uj 



Cg LU + iKS 2 /"f 



2nu 



mVl, 



2ttLU 



2 mVl 



(A9) 



mVl 



I W LU 



w 2 \ru n H luj + iks 2 



(A2) 



Consider first the case 51 = 0. Then (A9) is a quadratic 



equation in u> . Instability exists if there are wavenumbers 
for which lu 2 < 0. The sufficient and necessary condition 
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for instability is that the constant term in the quadratic 
is positive. The cases m = and m ^ have to be consid- 
ered separately. For both cases, the range of instability is 
maximized in the limit in which the vertical wavenumber 
n — > oo, and the instability conditions are 



p > 1 (m = 0), 



m 

P> T -1 



(m ^ 0). (A10) 



These are the conditions found by Tayler (1973) and 
Goossens et al. (1981), if we take into account that their 
results are valid on the entire sphere, whereas we have 
considered only the situation near the pole. 

For an azimuthal field resulting from the winding up of 
a a radial field component by differential rotation, B ~ w 
near the pole, so that the m = mode is marginally stable. 
The higher m's require a steeper field gradient than m = 1, 
hence m = 1 is the most unstable mode, at least in some 
region around the pole. From now on I ignore the m = 
mode. 

The case of arbitrary rotation rate is a bit complicated, 
so I consider only the limiting cases O = and Q/u>a — * 
oo. For O = 0, the dispersion relation can be made real 
by the substitution 

lo = ia. (All) 
Let K = ks 2 /7 and H = rjs 2 . Multiplying (AE) by m 2 (a + 



K)(a + H), one gets a fifth degree polynomial equation in 



-lo\( P - l)o-(<7 + H){a + K) + —W(cj + H)(a + K)+ 

I 2 

uj 2 A (a + K) + -^N 2 (a + H)}[a{a + H) + uj\) 



-uj 4 A (a + K)(p + 1) = 0. (O = 0) 



(A12) 



The system can in principle have both monotonic and os- 
cillatory (overstable) instabilities, and Eq. ( A12 ) would 
have to be checked fo r bot h possibilities. For systems with 
real coefficients like (A12), experience from double diffu- 
sive systems shows that overstability results if the destabi- 
lizing agent diffuses faster than the stabilizing agent, while 
monotonic instability results if the destabilizing agent has 
the lower diffusivity. In the present case the destabilizing 
agent is the magnetic field, which has a lower diffusivity 
than the stabilizing thermal stratification, so we expect 
monotonic instability. Marginal stability then corresponds 
to a = 0, and the condition for instability is that the con- 
stant term in the polynomial be negative. This yields 



m 



1 



I 2 jN 2 T) 
n 2 oj 2 . K 



(A13) 



For a smooth field gradient, p ~ 0(1), instability is pos- 
sible only if the last term does not exceed order unity: 



> I 



,N 2 rj 



(A14) 



This can be achieved, for arbitrarily low field strength, by 
taking the radial wavenumber n sufficiently large. There is 
no critical field strength for instability. Magnetic diffusion 
does have an effect on the growth rates, however. Closer 



inspection of Eq. (A12) shows that the maximum growth 



rate as a function of n is of the order of the adiabatic rate 
wa only if a critical field strength ujac is exceeded, given 
by: 



2 

w Ac 2~r iV ) 



(A15) 



where we have taken I ~ 1/r as the lowest possi- 
ble horizontal wavenumber. If this satisfied, the radial 
wavenumber at maximum growth rate is of the order 
n ~ (loa/t]) 1 ^ 2 . If the field strength is less than given by 



( /VI 5), the maximum growth rate is reduced: 



CT~uj a (wa » WAc ), <7~w|— (w A < W A c)-(Al6) 

Next consider the opposite limiting case, O ^> uj a . In 
this limit, one finds that the frequency scales as ui ~ w A /0. 
Writing 



lu — ao; A /0, h= T ^Y, k = 



2 ' 



(A17) 



and neglecting higher order terms in lua/Q, relation jAq ) 
reduces to 

m(p — l)(a + in 2 h)(a + in 2 k)+ 

m2 r -2, i2n2 / • 2, m 
——\a + in AH — s — ^-{a + in h)\ — 
2 n 2 uj A 

[2m(a+in 2 h)+p+l](a + in 2 k)[m(a + in 2 h) + l] = 0.(A18) 
Consider first the adiabatic case, h = k = 0. The dis- 



persion relation ( A18 ) is then quadratic in a and one finds 
that the necessary and sufficient condition for instability 



m 



(m ^ 0, Q > oj a ). 



(A19) 



This condition is significantly more restrictive than nonro- 
tating condition ( AlC| ). For a field B ~ w, such as would 
result from differential rotation near the pole, the condi- 
tion predicts stability. Since the fields we envisage are just 
of this type, one would conclude stabilily, at least in the 
interesting region near the pole that critical for driving 
the instability. This was noted by Pitts & Tayler (1986), 
who also found that a sufficiently large gradient p would 
again be unstable in the rapidly rotating case. It turns 
out, as shown next, that the instability condition for the 
rapidly rotating case is relaxed again when the effects of 
thermal and magnetic diffusion are taken into account. 
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Returning to Eq. (Al£) one easily verifies that there is 
no direct instability (a imaginary), since the coefficients 
are complex due to the combined effects of rotation and 
diffusion. Thus we have to check the stability of oscilla- 
tory modes. The stability boundary for marginally stable 
oscillations is found by requiring the dispersion relation to 
have a solution with a real. The case with both diffusiv- 
ities present is rather complicated, so I specialize further 
to the cases k = and k/t] — ► oo. The first limit is ap- 
propriate for cases where the stratification iV 2 is due to 
a composition gradient instead of the thermal gradient, 
because the ions making up such a gradient diffuse only 
slowly. The latter case applies when composition gradients 
can be ignored. 

For the case k — the real and imaginary parts of 
( [All ) yield 



Re 



Im 



-2(ma+l) 2 + 2+^-(p+l)- 
"'C r + 2m 2 n 4 h 2 = 0, 



2 n 2 ui 



-2m 2 a 2 — 2ma 



4 TL^iO 

(« = o, n > ua) 



T =0, 



Eliminating a between the real and imaginary parts yields 



(l + ia!44)l/2_l]2 + m!_ (p+1)+2TO 2 n 4 ft 2 = Q.(A20) 



The terms involving the wavenumber n are both positive. 
The last one increases monotonically with wave number, 
the first with the inverse of the wavenumber. Instabil- 
ity, if it exists, is therefore restricted to a finite range 
in wavenumbers. At the high wavenumber end the insta- 
bility is cut short by magnetic diffusion (last term), at 
low wavenumbers by the stable stratification (first term). 
I nsta bility exists if there are wavenumbers for which Eq. 
(A2C) has solutions. Necessary and sufficient for this to be 
the case is that the sum of the first and last terms be less 
than p + 1 — m 2 /2. This yields rather complicated expres- 
sions. Instead of the exact conditions, a sufficient condition 
for instability which is also sufficiently accurate as a neces- 
sary condition is found by noting that [(1+IE 2 ) 1 / 2 — 1] < X 2 , 
hence instability is guaranteed if there are n's for which 
both 



l 2 N 2 



<l(p + l_!£) 



(A21) 



where 

a = p + 1 - m 2 j2. 



(A24) 



With m = 1, I ~ 1/r, a nd a ~ 1, this is equivalent to the 
heuristic result of Sect. 5/2 (eq [l7|) . 

In th e opp osite limit n 3> rj, the real and imaginary 
parts of ( |Al|) yield 



2 - 2{ma + l) 2 + f 2 + \ - (p + 1) + g 2 = 0, 



-2g 2 /(f 2 + 2g 2 ), 



where 

f = 



, _ m 2 l 2 N 2 r] 
2 n 2 uj\ k ' 



g 2 = 2m 2 n 4 h 2 



The second of these implies that 
2 - 2(ma + l) 2 > 0. 



(A25) 



(A26) 



(A27) 



Necessary for instability is again that a — p+l—m 2 /2 > 0. 
Supposing that this is satisfied, a sufficient for instability 
is that there exist wavenumbers for which 

2-2(ma+l) 2 < -a, and f 2 < -a, and g 2 < -a.(A28) 

o o o 



This is the case if both p > 



1 and 



c4;>^ 2 iv 2 V^ 



(r) <C k, ft > ll>a)- 



(A29) 



This is a slightly imprecise condition, since I have not de- 
termined the exact numerical factor in front of the RHS. 
With to = 1, I ~ 1/r, and a ~ 1, this condition is equiva- 
lent to the heuristic result of Sect. 5.2 (eq ^9|). 



and 



2m 2 n 4 h 2 < |(p+l- 



(A22) 



This is possible only if the adiabatic instability condition 
p+ 1 — m 2 /2 > is satisfied, hence magnetic diffusion re- 
stricts the range of instability (as opposed to other double 
diffusive systems, where diffusion can be destabilizing). 

2 

Both conditions can be satisfied if both p > ^ 1 and 



<4 > 



a 3/2 



l 2 N 2 rjn, (k = 0, ft > lu a ) 



(A23) 
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